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I. INTRODUCTION 


Parametrically driven oscillators are of fundamental in¬ 
terest as a platform for studying quantum phenomena 
away from thermal equilibrium in well-characterized sys¬ 
tems. Examples of such phenomena include noise squeez¬ 
ing, low-noise quantum amplification, photon generation, 
quantum measurements, and coherent state preparation 
in a dissipative environment M- Classically, paramet¬ 
rically excited vibrational states are states with a bro¬ 
ken discrete time-translation symmetry, as their period is 
twice the modulation period [91] . Excitation of period two 
vibrations with the increasing driving strength is some¬ 
what reminiscent of the mean-field ferromagnetic phase 
transition with lowering temperature. 

A distinctive feature of the parametric oscillator is 
the possibility to have a stable state with unbroken 
time-translation symmetry along with the stable period- 
two vibrational states (more precisely, such states are 
called asymptotically stable). In the symmetric state, 
vibrations are either not excited or there may be small- 
amplitude vibrations with the same period as the driving. 
In the parameter space there is a critical point where this 
state and the period-two states merge. 

In this paper we study quantum fluctuations near the 
critical point. They are strong and slow. We find the 
long-time behavior of their correlation functions and the 
scaling of the correlation time with h and temperature. 
As the system moves deeper into the range of coexisting 
period-two states, in the neglect of fluctuations it would 
localize in one of them or in the symmetric state, if it is 
also stable. A major manifestation of quantum fluctu¬ 
ations is switching between the states. We identify the 
switching mechanism and find the switching rates. We 
study how the rates scale with the driving amplitude and 
how there occurs a crossover to the previously explored 
scaling far from the critical point. 

We focus on quantum effects for weakly damped oscil¬ 
lators, where the decay rate T is small compared to the 


oscillator eigenfrequency uiq- This condition is usually 
met in the experiments with superconducting microcav¬ 
ities and with Josephson-junction based systems used in 
quantum information. For weakly damped oscillators, 
nonlinearity becomes substantial already for compara¬ 
tively small vibration amplitudes, once the nonlinearity- 
induced shift of the vibration frequency becomes compa¬ 
rable with r. It is the nonlinearity that leads to multi¬ 
stability of forced vibrations. Even though the nonlinear 
effects are strong, they occur in the range of compara¬ 
tively small vibration amplitudes, and classically, the os¬ 
cillator motion is almost sinusoidal vibrations with slowly 
varying amplitude and phase. 



FIG. 1. The relation between the parameters of the driving 
field where the number of stable vibrational states changes. 
Parameters and / p characterize the scaled detuning of the 
field frequency from resonance and the scaled field amplitude, 
respectively, see Eq. ©. For f p < 1 or fi p < p,Bi the oscillator 
has no period-two vibrational states. In the region / p > 1 
and fip > pbi there are two stable period-two states. For 
Mp > 9 B 2 the state with no period-two vibrations is stable, 
too. On the long-dashed line the stable and unstable period- 
two states coalesce. The point fi p = 0, / p = 1 where the three 
bifurcation lines merge is the critical point. The dot-dashed 
line shows where the three stable states are equally occupied 
in the stationary regime. 

Classically, the dynamics of a driven underdamped os¬ 
cillator is determined by two parameters, the scaled fre- 
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quency and amplitude of the driving [ 9 ]. The analysis 
of the dynamics is simplified near bifurcation parameter 
values, where the number of the stable vibrational states 
changes. In this region one of the motions becomes slow, 
an analog of the “soft mode” 0- The bifurcation re¬ 
lation between the oscillator parameters in the case of 
close-to-resonant parametric driving is shown in Fig. |T] 

The possibility to describe the dynamics by a single 
dynamical variable also significantly simplifies quantum 
analysis [Tlj| . Such variable commutes with itself at dif¬ 
ferent times, and therefore its dynamics is essentially 
classical. However, the fluctuations are quantum, be¬ 
cause the noise that causes them has quantum origin and 
its intensity is oc h for low temperature. As we show, 
this picture applies near the critical point. We note that, 
in terms of the Floquet states of the modulated oscilla¬ 
tor, this case corresponds to the distance between the 
quasienergy levels being much smaller than their width 
due to dissipation. 

II. THE MODEL 

We consider a parametrically driven oscillator with the 
Duffing (Kerr) nonlinearity. The Hamiltonian of the os¬ 
cillator reads 

H 0 = + ^q 2 F cos uj F t, ( 1 ) 

where q and p are the oscillator coordinate and momen¬ 
tum, the mass is set equal to one, and 7 is the anhar- 
monicity parameter. Parameter F gives the amplitude 
of the driving, whereas ui F is the driving frequency. 

We assume the driving to be resonant and not too 
strong, 

5 u)=^uj f -u) 0; |7l <9 2 )<Wq. (2) 

It is convenient to change to the rotating frame 
using the standard canonical transformation U(t) = 
exp (— 10 s au>Ft/ 2 ), where a 1 and a are the raising and 
lowering operators of the oscillator. We introduce slowly 
varying in time dimensionless coordinate Q and momen¬ 
tum P, using as a scaling factor the characteristic am¬ 
plitude of forced vibrations C p = \2F C / 3 -y | 1//2 [parameter 
F c is defined below in Eq. ©], 

U\t)qU(t) = C P (Q cosip + Psinip), 

U^(t)pU(t) = — -C p ujf(Q sirup — P cosp) ( 3 ) 

with ip = ^u>pt + j7r. The commutation relation between 
P and Q is 

[Hi Q] = — -V = 3\'y\fc/ (4) 

Parameter A p oc h plays the role of the Planck constant 
in the quantum dynamics in the rotating frame. It is 
determined by the oscillator nonlinearity, A p oc 7, For 


characteristic |Q|,|P| < 1 , where (q 2 ) < C p , the last 
inequality in Eq. © implies F <C To simplify the 
analysis near the critical point, in Eqs. © and we 
choose the variables in the form that slightly differs from 
that in Ref. nn 

In the range © the oscillator dynamics can be ana¬ 
lyzed in the rotating wave approximation (RWA). The 
Hamiltonian in the rotating frame is U^HoU — ihUW = 
(P c 2 /67)<? p . Operator g p = g p {Q,P) is independent of 
time and has the form 

g P (Q,P) = \ (p 2 + Q 2 ) 2 - \v P (P 2 + Q 2 ) 

+ ^/pSgn7(QP + PQ), 

u) F (ujf — 2w 0 ) , , _s 

Pp = -F-sgny, / p = F/F c . ( 5 ) 

F C 

The eigenvalues of g p (multiplied by F 2 / 67) give the os¬ 
cillator quasienergies, 

In Eq. ©, parameter / p is the scaled driving amplitude 
P; in what follows we assume / p > 0 . Parameter g p 
gives the detuning of the driving frequency from twice the 
oscillator eigenfrequency. The scaling factor in both f p 
and /i p is the critical amplitude F c needed for parametric 
excitation of the oscillator. 

Relaxation of the oscillator results from the coupling 
to a thermal bath. The coupling leads to oscillator decay 
due to scattering by the bath excitations. For a weak cou¬ 
pling, one of the most important scattering mechanisms 
is scattering with energy transfer « Huio in an elementary 
event Il 2 - [l 4 | . It comes from the coupling, which is lin¬ 
ear in the oscillator coordinate and/or momentum. In a 
phenomenological description of the oscillator dynamics 
such scattering corresponds to a friction force — 2 Tq. The 
friction coefficient T is simply expressed in terms of the 
appropriate correlator of the thermal bath variables. 

It is important for what follows that relaxation comes 
along with a quantum noise. The equations of motion for 
slow variables Q, P in the rotating frame in dimensionless 
time r = Ft read 

Q = -AAp 1 ^, <? p ]sgn7 - Q + /q(t), 
P=-i\- 1 [P,g p ]sgn^-P + fp(T). (6) 

Here, A = dA/dr. We have set 

F c = 2 Toj f , ( 7 ) 

in which case, indeed, F c is the threshold value of the 
driving amplitude, as will be seen below. 

In Eq. ©, f q p are quantum noise operators. The 
noise is ^-correlated in slow time, cf. [Tli] , 

(/qM/q^')) b = (/p(r)/p(r , ))b = 2 V6(t - r'), 

([/q(t). /p( T ')])b = 2 iX p S(r - r'), 

V = A p (h + 1/2), n = [exp^hujo/ksT) — l] -1 . (8) 
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Here, n is the oscillator Planck number and (-)b means 
averaging over the bath state. Parameter V plays the 
role of the effective temperature of the quantum noise, 
T> oc h for hut o fc B T. The commutation relation in 

Eq. J8]) guarantees that the commutator [Q, P] = i\ p 
does not change in time. The noise correlators are un¬ 
derstood in the Stratonovich sense []1|; in particular, 

([fQ( T ),P( T )])b = ([<9(r),/p(r)])b =*A P . 

III. SLOW QUANTUM DYNAMICS NEAR THE 
CRITICAL POINT 

The stable vibrational states of the oscillator are given 
by the stable stationary solutions of Eq. (JG]) without 
noise. The bifurcation values of the parameters / p ,/i p , 
where the number of the stable states changes, lie on the 
lines p p = ±/i B (/ p ) and (p p > 0, / p = 1) on the (/ p , p p )- 
plane, 

Mb = (/ p 2 - 1) 1/2 - (9) 

The bifurcation lines are shown in Fig. [I] For p p < psi = 
~Pb and for f p < 1 the only stable state of the oscillator 
is Q = P = 0. The amplitude of period-two vibrations in 
this state is zero, and below we call it the zero-amplitude 
state. At /i p = p B \ this state becomes unstable. For 
larger fi p or f p the oscillator has two stable states, which 
correspond to period-two vibrations in the lab frame with 
the reduced squared amplitude Q 2 + P 2 = /i p + (f 2 — 
l) 1 / 2 and with the phases differing by n. For /i p = /i B2 
the state Q = P = 0 becomes stable and there emerge 
two unstable states, which correspond to unstable period- 
two vibrational states in the lab frame [9]. The values of 
the control parameter /i B i and /x B2 correspond to the 
supercritical and subcritical pitchfork bifurcations Ed- 
On the line f p = 1, /i p > 0 the stable and unstable states, 
which correspond to period-two vibrations, coalesce. 

At the critical point (/r p = 0, / p = 1) all three bifur¬ 
cation lines merge. Such merging is robust, it does not 
disappear if the model is slightly changed, although the 
position of the critical point can change. The associated 
strong singularity of the dynamics in the absence of noise 
leads to comparatively strong fluctuation effects. 

Near the critical point, for small Q, P the coordinate Q 
varies in time much slower than P. This is seen from the 
linearized equations ([6]) , which in the absence of noise 
take the form Q = (/ p — 1 )Q — p p Psgn~/, P = — (/ p + 

1 )P + MpQsgny. 

Because Q(t) is slow, function P(t) adiabatically fol¬ 
lows Q(t). For the dimensionless time r>l, function 
P[t) can be expressed in terms of Q(t) by setting P = 0. 
Substituting the result into the full nonlinear equation for 
Q , one obtains 

Q w -d Q U(Q) + f Q (T), (10) 

U{Q) = -Q 2 [m p - (f P - 1)] + jQ 4 ^-m p + 3*3^ 


[U ( Q ) has an extra factor 2/(/ p +1) which is set equal to 
1, to the leading order in / p — 1]. We note that, because 
the relaxation time of P(r) is small compared to that 
of Q(t), fluctuations of P(t) due to the quantum noise 
/ B (r), are small compared to fluctuations of Q(t). This 
is the squeezing effect. 

For psi < p P < MB2, the effective potential U(Q) has 
a local maximum at Q = 0 and two symmetric minima 
at Q q 2 = ±(/r B + Mp ) 1 ^ 2 ■ The minima correspond to the 
stable period-two vibrations in the lab frame. 

For n p > HB2, f p > 1, function U(Q) still have the 
minima at Q i >2 , but now it has a local minimum at Q = 0 
and two additional local maxima at Q = ±(/x p — ^lb) 1 2 ■ 
The minimum at Q = 0 corresponds to the stable zero- 
amplitude state. The local maxima correspond to unsta¬ 
ble period-two vibrations. 

The behavior of the oscillator near the onset of period- 
two vibrations at /x B i has similarities with the mean-field 
picture of the critical behavior at the second-order phase 
transition. The time-translation symmetry of the stable 
stationary state at Q = 0 is spontaneously broken and 
there emerge two vibrational states with equal ampli¬ 
tudes and opposite phases and with period 47r/wp. Such 
behavior is usually described by an effective potential 
which is quartic in the coordinate Q. 

A key observation is that, near the critical point f p = 
1,/Xp = 0, it is necessary to keep the sixth-order term 
in U(Q). This follows from Eq. (flOl) . The system be¬ 
comes softer than at the bifurcation point far away from 
the critical point. Such form of the potential reminds 
the Landau free energy for a system that can undergo a 
first-order phase transition, and indeed, an analog of this 
transition can occur in the driven oscillator, see below. 


A. The quantum Fokker-Planck equation 

Equation m has the form of the Langevin equation 
of a classical particle with one dynamical variable (half¬ 
degree of freedom) driven by a <5-correlated noise. Indeed, 
/q(t) and Q(t) commute, and therefore [Q(t),Q(t')\ = 
0. However, the noise has quantum origin, its correla¬ 
tor flSJ) explicitly contains H. The slow quantum fluc¬ 
tuations can be alternatively described by the Fokker- 
Planck equation for the probability density p(Q, r). From 
Eq. ifTOl) . it reads 

d T p = d Q (pd Q u) + vd 2 Q p. (ii) 

The stationary distribution of the oscillator near the 
critical point has the Boltzmann form 

p st (Q) = Z- 1 exp[-U(Q)/V\, (12) 

where Z = J dQ exp[— U(Q)/T>] is the effective partition 
function. The evolution of this distribution with varying 
f p is shown in Fig. [2] 
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FIG. 2. Evolution of the stationary distribution p B t(Q) with 
varying amplitude of the driving field for fi p = 0. For f p < 1 
the distribution has a single broad maximum at Q = 0. With 
increasing / p there are developed two sharp maxima at the 
positions that correspond to the period-two vibrational states. 

The full stationary Wigner distribution is a product of 
Pst(Q) and a Gaussian distribution over P, with the po¬ 
sition of the maximum over P being Q-dependent. Func¬ 
tion Pst (Q) can be thought of as the integral over P of 
this full distribution. The distribution m is maximal 
at the minima of U(Q), i.e., at the values of Q that cor¬ 
respond either to the stable states of the oscillator. Near 
the critical point the distribution is broad, which corre¬ 
sponds to large fluctuations in this range. 


IV. SCALING OF THE RATES OF 
QUANTUM-FLUCTUATION INDUCED 
SWITCHING 

Once the depth of the minima of U ( Q ) becomes large 
compared to A p (2h +1), the dynamics of the oscillator 
is characterized by two very different time scales. One is 
the relaxation time l/U"(Q) for Q near the stable states. 
The other is a much slower time over which the oscilla¬ 
tor can switch between the minima due to comparatively 
rare large fluctuations. The switching rate Ws W is given 
by the Kramers theory |lj| of thermally activated switch¬ 
ing over a potential barrier. This theory immediately ex¬ 
tends in the present case to quantum fluctuations, even 
though there is, of course, no thermal activation for low 
temperatures. In dimensionless time t = Tt 

W sw = Q sw exp(—P^/Ap), Ra = At//(h + 1/2), 

aw = [|C/"(Qs)|C/"(Oa)] 1/2 /27r, 

AU = U(Q s )~U(Q a ). (13) 

Here, Q a is the position of the stable state (attractor) 
from which the oscillator switches and Qs is the posi¬ 
tion of the adjacent local potential maximum [a saddle 
point on the ( Q , P)-plane] over which the switching oc¬ 
curs. Function Ra is the quantum activation energy for 
the switching, it is proportional to the height of the po¬ 
tential barrier overcome in the switching. 


Remarkably, Eq. m for the switching rate has the 
same structure as the expression for the rate of tunnel¬ 
ing from a potential well of U(Q) in the sense that the 
rate is exponential in 1/A P oc 1/S. However, switching 
occurs not via tunneling under the barrier of U(Q) but 
by going over the barrier. Tunneling requires the fast 
variable P to be involved. Therefore the tunneling expo¬ 
nent is parametrically larger than Ra and the tunneling 
rate is exponentially smaller E3- 

The states between which the switching occurs and 
the rate FF SW depend on the parameter region on the 
bifurcation diagram in Fig. [T] We start with the region 
PBi(fp) < Pp < PB 2 (fp )• Here, the oscillator is bistable. 
It can switch from one period-two state to the other. By 
symmetry, the switching rates for the both states are the 
same. In Eq. m Q a = ±{pb + Sp) 1 ^ 2 ; depending on 
the initially occupied state. The switching occurs over 
the unstable state Qs = 0. 

The dependence of the activation energy Ra for switch¬ 
ing between the period-two states on the parameters is 
shown in Fig. [3] For exact resonance, p p = 0, it fol¬ 
lows from Eq. (fTOl) that Ra = (/ p — l) 3 / 2 /3(2h + 1). 
In this case Ra displays a power law dependence on the 
distance from the critical field amplitude f p — 1. The 
scaling exponent is 3/2. For nonzero p p , the behavior be¬ 
comes more complicated and is not described by a simple 
power law. Close to the bifurcation line pBi(fp), func¬ 
tion Ra scales as a power of the distance from p.Bi(fp), 
Ra ~ Pb{p p — Msi) 2 /2(2h + 1), cf. Ref. [U We note 
that the scaling of the switching rates near a bifurcation 
point of a different type has been observed for resonantly 
driven Josephson junction based quantum oscillators [18| ■ 
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FIG. 3. The quantum activation energy Ra = (n + 1/2 )Ra 
for switching between the parametrically excited period-two 
vibrational states with opposite phases. In the shown range 
of the scaled driving field amplitude 1.04 < f p = F/F c < 1.2 
and the scaled frequency detuning |/r p | < /rs(/ p ), the period- 
two vibrations are the only stable states of the oscillator. 

A different situation occurs in the range p p > pB 2 (fp) 
and f p > 1. Here, far from the bifurcation lines in Fig. [T] 
the potential U ( Q) has three local minima separated by 
two local maxima. Interstate switching occurs via a tran¬ 
sition from a minimum of U ( Q ) over the adjacent maxi¬ 
mum. This means that from the stable period-two states 
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the oscillator can switch to the zero-amplitude state and 
vice versa, but direct switching between the period-two 
states is exponentially unlikely. 

The switching activation energies in the region of 
tristability are shown in Fig. 0] Rai and Rao refer to 
switching from a period-two state to the zero-amplitude 
state and from the zero-amplitude state to one of the 
period-two states, respectively. To find Rai, one sets in 
Eq. (fldl) Q a = ±(/r p + yiis) 1 / 2 , whereas to find Rao one 
sets Q a = 0; in the both cases, Qs = ±(/i p — P-b ) 1 ^ 2 ■ In¬ 
terestingly, from Eq. (in, the activation energy Rai dis¬ 
plays simple scaling behavior, Rai = (/ p — l) 3 / 2 /3(n + 
1/2) independent of p p . The activation energy of switch¬ 
ing from the zero-amplitude state displays a more com¬ 
plicated behavior. It decreases with the increasing field 
amplitude. Near the bifurcation line p p = pB(fp) we 
have Rao ~ Uniu-n — uft?) 2 /2(2n+ 1), in agreement with 
the previous work [ll| . 
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FIG. 4. The quantum activation energies -Rai = (n+l/2)RAi 
and Rao = (n+ 1/2) Rao for switching from a state of period- 
two vibrations to the zero-amplitude state and from the zero- 
amplitude state to one of the states of period-two vibrations, 
respectively. The curves refer to the region of tristability, 
where the scaled driving field amplitude f p > 1 and the scaled 
frequency detuning p p > PB(f P )- 


In the stationary regime, for Rai > Rao it is ex¬ 
ponentially more probable for the system to be in the 
state of period-two vibrations. On the other hand, for 
Rao > Rai, preferentially occupied is the zero-amplitude 
state, where the period-two vibrations are not excited. 
The interrelation between the field amplitude and fre¬ 
quency where Rai = Rao corresponds to an analog 
of a first-order phase transition: here the populations 
of the different stable states of the oscillator are close 
to each other. From Eq. m, near the critical point 
the “phase-transition” value of f p is given by expression 
(fp — l) 1 / 2 = /ip/2. It is shown by the dot-dashed line 
in Fig. □ 


V. QUANTUM FLUCTUATIONS IN THE 

IMMEDIATE VICINITY OF THE CRITICAL 
POINT 

Close to the critical point, where |/ 2 —1| < D 2 / 3 , |/z p | < 
I? 1 / 3 , the stationary distribution (1121) is almost flat, as 
seen from Fig. [2] This indicates anomalously large quan¬ 
tum fluctuations of the oscillator variable Q (whereas the 
fluctuations of the other variable, P, are squeezed). The 
fluctuational mean-square displacement is ( Q 2 ) ~ D 1 / 3 , 
whereas far from the critical region the mean-square dis¬ 
placement about a stable state is ~ V (V oc h is the small 
parameter of the theory, D C 1). In the critical region 
the exponents Ra/A p are no longer large and the con¬ 
cept of switching rates becomes ill-defined. One cannot 
separate switching from other fluctuations. 

A convenient characteristic of the fluctuation dynam¬ 
ics is the long-time decay of the correlation functions of 
the oscillator. This decay is characterized by the low¬ 
est nonzero eigenvalue iq of the Fokker-Planck equa¬ 
tion m • It can be found in a standard way fTbt | by 
reducing Eq. m to a Schrodinger-type equation for 
p(Q) = exp\U(Q)/2V\p(Q). The quantum noise inten¬ 
sity V can be scaled out of this equation by changing to 
the scaled time P 2 / 3 r, scaled coordinate r D~ 1 / 6 Q, and 
scaled parameters V~ 1 / 3 p p and P“ 2 ^ 3 (/ p — 1). 

The dependence of the decrement V\ on the frequency 
detuning parameter p p oc (wp — 2wo)/P c for several values 
of the scaled driving field amplitude / p = F/F c is shown 
in Fig. [5] It is seen that, for a given p p , the dynamics is 
slowed down, i.e., zq decreases with the increasing driving 
amplitude. This is to be expected and is in agreement 
with Fig. [2] For / 2 < 1 or (/ 2 - l) 1 / 2 + p p < 0, the 
oscillator has one stable state, Q = 0. As / 2 increases, 
the minimum of U(Q) at Q = 0 becomes more and more 
shallow and the relaxation rate, which is related to the 
curvature of U(Q ), goes down. As / 2 further increases, 
the system goes into the regime of coexisting period-two 
states (via the range of coexisting three stable states, 
for p p > 0, cf. Fig. [[]). For large (/ 2 — 1)/D 2 / 3 , zq is 
determined by the exponentially small rate of interstate 
switching that rapidly decreases with increasing / p , see 
Figs. [3] and 01 

The dependence of zq on the frequency detuning in 
Fig. [5]is profoundly nonmonotonic. For / 2 —1 < 0, where 
the system has only one stable state, Q = 0, this behav¬ 
ior is a consequence of the nonmonotonic dependence of 
the curvature of U"( 0) on p p , see Eq. (ITOl) . Close to the 
critical point the nonparabolicity of U ( Q) is substantial. 
This is also seen from Fig. (2j which shows that the sta¬ 
tionary probability distribution becomes profoundly non- 
Gaussian near the critical point. Therefore iq is different 
from U"( 0). The minimum of v\ is reached for p, p > 0. 

In the range p p > ps l = ~(f p — l) 1 ^ 2 the oscillator 
has two or three stable states in the neglect of fluctua¬ 
tions. For large (p p — Pbi)/F > 1 / 3 , the decrement iq is 
determined by the rate of interstate switching. This ex¬ 
plains the decrease of zq with increasing p p in Fig. [5] for 
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FIG. 5. The smallest nonzero eigenvalue v\ of the kinetic 
equation m, which determines the long-time decay of the 
correlation functions very close to the critical point = 
0, f p = 1. The curves from top to down correspond to 
(/p 1)/X? 2 / 3 = -4, 2, 0,2,4,6. 

/p > 1 and for not too large n P /V 

To understand the behavior of v\ for larger fi p in the 
range of tristability of the oscillator, we consider the bal¬ 
ance equations for the populations of the two period-two 
states 101,2 and the zero-amplitude state wq, 

w 0 = -2W 0 iw 0 + Wio(wt + w 2 ), 

Wi = Woitoo - W w wi (i = 1,2), (14) 

where Woi = W / o 2 oc exp(— Rao/\) is the rate of switch¬ 
ing from the zero-amplitude state to one of the period- 
two states and Wio = W 2 q oc exp(— Rai/X p ) is the rate of 
switching from a period-two state to the zero-amplitude 
state. 

From Eq. m v\ = W\q in the region of tristabil¬ 
ity where the slow dynamics is determined by inter-state 
switching. The exponent Rai/\ in Eq. (THU) for the rate 
Wio is independent of /U p , cf. Fig. |4j However, the pref¬ 
actor O sw oc (/ 2 - l) 1/2 [/ip - (fp - 1)] 1/2 is increasing 
with increasing /r p . Therefore, somewhat unexpectedly, 
v\ starts slowly increasing with /i p in the tristability re¬ 
gion. 

VI. CONCLUSIONS 

The results of this paper refer to quantum fluctuations 
near the threshold of parametric excitation of a nonlinear 
oscillator. The threshold corresponds to the frequency of 
the driving field u>f equal to twice the oscillator eigen- 
frequency ui o and to the field amplitude F equal to the 
critical value F c proportional to the oscillator decay rate. 
For small | F — F c |, |wf — 2u>o\, the oscillator dynamics is 
mapped onto the dynamics of an overdamped Brownian 
particle driven by quantum noise and confined in a sym¬ 


metric potential well of the sixth order in the particle 
displacement. 

In the immediate vicinity of the threshold quantum 
fluctuations are anomalously large and slow. The root- 
mean-square fluctuations of the amplitude scale with h 
and with temperature as [h(2h + l)] 1 / 6 . The correlation 
time of the fluctuations scales as [H(2h + l)]^ 2 / 3 . We 
have found the dependence of the correlation time on 
the frequency detuning uip — 2 ujq and the field amplitude. 
Unexpectedly, the dependence on the frequency detuning 
turned out to be nonmonotonic. 

For F > F Cl in the appropriate frequency range the 
oscillator can have coexisting states of period-two vibra¬ 
tions with the interstate distance in phase space sig¬ 
nificantly exceeding the mean-square-root fluctuations 
about the states, which are ~ [h(2fi + l)] 1 / 2 . Here, rare 
large quantum fluctuations lead to switching between 
the states. The switching mechanism is quantum acti¬ 
vation. In switching the oscillator goes over an effec¬ 
tive quasienergy barrier rather than tunneling beneath 
it. However, the logarithm of the switching rate is oc 1/h 
for low temperatures, as in the case of tunneling. 

We found simple explicit expressions for the switch¬ 
ing rates both in the regime where the period-two states 
are the only stable states of the oscillator and where the 
state with no period-two vibrations is also stable. In 
the region of bistability, at exact resonance the switch¬ 
ing activation energy scales with the field amplitude as 
(F — F c ) 3 / 2 . This scaling applies also to the switching 
from period-two states in the regime of tristability for 
arbitrary detuning. In contrast, the dependence of the 
activation energy on F — F c in the regime of bistabil¬ 
ity away from exact resonance, and the dependence on 
F—F c of the activation energy of switching from the zero- 
amplitude state do not display simple scaling, generally. 
In the region of tristability, depending on the field param¬ 
eters, either the period-two states or the zero-amplitude 
state are predominantly occupied. The state occupations 
are close in a narrow parameter range, the behavior that 
reminds the first-order phase transition. 
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